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ABSTRACT 

Recent advances in the understanding of the properties of supernova remnant 
shocks have been precipitated by the Chandra and XMM X-ray Observatories, 
and the HESS Atmospheric Cerenkov Telescope in the TeV band. A critical 
problem for this field is the understanding of the relative degree of dissipative 
heating/energization of electrons and ions in the shock layer. This impacts the 
interpretation of X-ray observations, and moreover influences the efficiency of 
injection into the acceleration process, which in turn feeds back into the thermal 
shock layer energetics and dynamics. This paper outlines the first stages of our 
exploration of the role of charge separation potentials in non-relativistic electron- 
ion shocks where the inertial gyro-scales are widely disparate, using results from 
a Monte Carlo simulation. Charge density spatial profiles were obtained in the 
linear regime, sampling the inertial scales for both ions and electrons, for differ- 
ent magnetic field obliquities. These were readily integrated to acquire electric 
field profiles in the absence of self-consistent, spatial readjustments between the 
electrons and the ions. It was found that while diffusion plays little role in mod- 
ulating the linear field structure in highly oblique and perpendicular shocks, in 
quasi-parallel shocks, where charge separations induced by gyrations are small, 
and shock-layer electric fields are predominantly generated on diffusive scales. 



1. Introduction 

The understanding of the character of shells and interiors of supernova remnants (SNRs) 
has been advanced considerably by groundbreaking observations by the Chandra X-ray Ob- 
servatory. These have been enabled by its spectral resolution coupled with its impressive 
angular resolution. Of particular interest to the shock acceleration and cosmic ray physics 
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communities is the observation of extremely narrow non-thermal X-ray spatial profiles in 
selected remnants (see Long et al. 2003 for SN1006; Vink & Laming 2003 for Cas A; and 
Ellison & Cassam-Chenai 2005, for theoretical modeling), which define strong brightness 
contrasts between the shell, and the outer, upstream zones. If the synchrotron mechanism 
is responsible for this non-thermal emission, the flux ratios from shock to upstream indicate 
strong magnetic field enhancement near the shock. These ratios considerably exceed values 
expected for hydrodynamic compression at the shocked shell, so proposals of magnetic field 
amphfication (e.g. Lucek & Bell 2000) in the upstream shock precursor have emerged. 

Another striking determination by Chandra concerns electron heating by ions in the 
shock layer. Dynamical inferences of proton temperatures in remnant shocks can be made 
using proper motion studies of changes in a remnant's angular size, or more direct spec- 
troscopic methods (e.g. Ghavamian et al. 2003). In the case of remnant IE 0101.2-7129, 
Hughes et al. (2000) used a combination of ROSAT and Chandra data spanning a decade 
to deduce an expansion speed. Electron temperatures Tg are determined using ion line di- 
agnostics (assuming the equilibration Te —Tp), via both the widths and relative strengths 
for different ionized species. From these two ingredients, Hughes et al. (2000) observed that 
deduced proton temperatures were considerably cooler, i.e. 3kTp/2 <^ mp(3'Ui^/4)^/2 , than 
would correspond to standard heating for a strong hydrodynamic shock with an upstream 
fiow speed of uix ■ The same inference was made by Decourchelle et al. (2000) for Kepler's 
remnant, and by Hwang et al. (2002) for Tycho's SNR. This property is naturally expected 
in the nonlinear shock acceleration scenario that is widely used in describing cosmic ray and 
relativistic electron generation in SNRs: the highest energy particles tap significant fractions 
of the total available energy, leading to a reduction in the gas temperatures. This nonlinear 
hydrodynamic modification has been widely discussed in the cosmic ray acceleration litera- 
ture (e.g. sec Jones & Ellison 1991; Berczhko & Ellison 1999, and references therein), and 
has been extensively applied to multiwavelcngth SNR spectral models (e.g. see Baring et al. 
1999; Berczhko ct al. 2002; Ellison & Cassam-Chenai 2005; Baring, EUison & Slane 2005). 

The extent of equilibration between electrons and ions in SNR shell shocks needs to be 
understood, and can potentially be investigated by laboratory plasma experiments. A critical 
ingredient is the electrostatic coupling between electrons and protons in the shock layer, 
which offers the potential for considerable heating of e~ , coupled with cooling of protons, 
setting mp{3uix/Af/2 > 3kTe/2 > me{3ui^/Af /2 with Te^Tp. Probing this coupling is 
the subject of this paper. Here we describe preliminary results from our program to explore 
electrostatic energy exchange between these two species in SNR shocks, using a Monte Carlo 
simulation of charged particle transport, their spatial distribution and associated electric 
field generation. The goal is to eventually obtain a simulation with self-consistent feedback 
between the charge separation potentials, and the Lorentz force they impose on the charges. 
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The research progress outlined here indicates that the role of diffusion in quasi-parallel shocks 
is very important, and can readily influence computed charge separation potentials. 

2. Shock Layer Electrostatics in Supernova Remnants 

Cross-shock electrostatic potentials arise in the shock layer because of the different 
masses of electrons and ions: upstream thermal ions gyrate on larger scales than do their 
electron counterparts when they transit downstream of the shock for the flrst time. In shocks 
where the fleld is oblique to the shock normal by some angle ©Bni upstream (and therefore 
a greater angle downstream), on average, protons will be located further downstream of the 
shock than electrons. This naturally establishes an electric field E, to which the plasma 
responds by accelerating electrons and slowing down ions to short out the induced E. A 
feedback loop ensues, mediated by fields and currents that vary spatially on the order of, 
or less than, the ion inertial scale, which is typically shorter than the ion gyroradius in 
astrophysical shocks such as those associated with supernova remnant shells. 

Particle-in-cell (PIC) simulations are a natural technique (e.g. sec Forslund & Friedberg 
1971 for an early implementation) for exploring signatures of such electrostatics in shock 
layers. These trace particle motion and fleld fluctuations, obtained as self-consistent solutions 
of the Newton-Lorentz and Maxwell's equations, in structured zones or cells in spatially- 
constrained boxes. Such simulations have been used recently to probe the Weibel instability 
at weakly-magnetized, perpendicular, relativistic pair plasma shocks (see Silva et al. 2003; 
Hededal et al. 2004; Nishikawa et al. 2005). They have also been used by Shimada & 
Hoshino (2000) to treat electrostatic instabilities at non-relativistic, quasi-parallel electron- 
ion shocks. Rich in their turbulence information, due to their intensive CPU requirements, 
such simulations have difficulty in modeling realistic nip/ rUe mass ratios, and fully exploring 
3D shock physics such as diffusive transport. Moreover, they cannot presently address the 
wide range of particle momenta and spatial/temporal scales encountered in the acceleration 
process; they often do not obtain time-asymptotic states for the particle distributions. 

Monte Carlo techniques provide an alternative method that can easily resolve electron 
and proton inertial scales, treat fully 3D transport and large dynamical ranges in spatial and 
momentum scales, all at modest computation cost. While they parameterize the effects of 
turbulence via diffusive mean free paths (e.g. see Jones & EUison 1991), they can accurately 
describe the microphysics of cross-shock electrostatic potentials. This simulational approach 
has been well-documented in the literature (e.g. Jones & Ellison 1991; Ellison, Jones & 
Baring 1996), with definitive contributions to the study of heliospheric shock systems, cosmic 
ray production, SNR applications and gamma-ray bursts. It models the convection and 
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diffusion of gyrating particles in spatially-structured flows and flelds, with transport back and 
forth across the shock effecting diffusive Fermi-type acceleration directly from the thermal 
population. The mean free path A is usually prescribed as some increasing function of 
particle momentum p or gyroradius Vg . Here we use this approach, with X cc p adopted as 
a broadly representative situation: see Baring et al. (1997) for a discussion of evidence from 
observations and plasma simulations in support of such a specialization. Here, A/r^ = 5 is 
chosen for illustrative purposes, to begin to investigate electrostatic influences on thermal and 
low-energy non-thermal particles in non-relativistic electron-ion shocks of arbitrary ©b„i . 
Clearly then, the diffusive scales for protons and electrons are disparate by their mass ratio. 

In the Monte Carlo simulation, the shock is defined magnetohydrodynamically, con- 
sisting of laminar, uniform flows and flelds upstream and downstream of a sharp, planar 
discontinuity. The magnetic flelds and flow velocities either side of the shock are related 
uniquely through the standard Rankine-Hugoniot solutions for energy and momentum flux 
conservation (e.g. see Boyd & Sanderson 1969). These solutions include essential elements of 
Maxwell's equations, such as the divergenceless nature of the B field. In the reference frames 
of the local upstream and downstream fiuids, the mean electric field is assumed to be zero, a 
consequence of very effective Debye screening, so that the only electric fields present in shock 
rest frames are u x B drift fields. The charged electrons and protons (more massive ions are 
omitted in this paper to simplify the identification of the principal effects) are treated as 
test particles, convecting into and through the shock, initially with the prescribed upstream 
fluid velocity Ui. This neutral beam is entirely thermal, and moreover is in equipartition, 
so that it has an input temperature T^. — Tp. The charges constantly diffuse in space to 
mimic collisions with magnetic turbulence that is putatively present in the shock environs, 
and in so doing, can be accelerated. These non-thermal particles form a minority of the total 
population, and provide only a minor contribution to the flelds illustrated in this paper. 

The charges transiting the shock distribute their downstream density in a manner that 
couples directly to their gyrational motion (e.g. see Baring 2006), and the local densities of 
electrons and protons can easily be tracked in the Monte Carlo technique by accumulating 
"detection" data at various distances from the shock. Monte Carlo simulation runs clearly 
exhibit non-zero charge excursions within a proton gyroradius of the shock, an effect similar 
to those found in PIC codes. For example, a cold, neutral e — p upstream beam develops 
an electron concentration near the shock in the downstream region, with protons distributed 
on their larger inertial scales. The resulting charge distributions p{x) depend on both the 
upstream fleld obliquity ©Bni , and also on the sonic Mach number M.^ ^ ui^/ y/^hTp/ (3mp) 
in situations where the upstream beam is warm. Due to the steady-state, planar nature of 
the simulation, these distributions depend only on the coordinate x along the shock normal. 
It is straightforward using Gauss' law for electrostatics, V • E = Airp^x) , to integrate the 
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x/(rjmpUjj.c/eBj) - Distance normal to shock, x x/(rjmpUjj^c/eBj) - Distance normal to shock, x 

Fig. 1. — Electrostatic profiles for an electron-proton plasma shock of speed uix = 9000 km/sec 
and sonic Mach number ^As = 10 , with upstream field obliquities 0Bni = 0° (green), 0Bni = 30° 
(blue) and ©Bni = 60° (red), as labelled. Left panel: scaled charge density distribution p{x) ; Right 
Panel: resulting ^^lineai^^ electric field profile E = Ex{x) computed by solving Gauss' Law (the 
©Bni = 0° displays Ex{x)/2). The profiles are exhibited on proton gyroscales Vg^p = TUpUixc/eBi , 
so that fluctuations on electron gyroscales are collapsed into the shock layer (dashed vertical line). 
The panels display both oscillations associated with proton gyrations in oblique cases, and diffusive 
upstream "precursors," which are most prominent when 0Bni = 0° (a parallel shock). 

charge distribution profile to obtain Ex{x) = —d^/dx. Eventually, such ^^linear'^ fields will 
then be used to compute the energy exchange between electrons and ions as they cross the 
non-monotonic charge separation potential $(x) . 

Linear determinations of electrostatic spatial profiles are shown in Figure 1 to illustrate 
the key features; these did not self-consistently include the acceleration of electrons and pro- 
tons in the produced E field. The left panel depicts large charge density fluctuations that 
trace ion gyration in the downstream magnetic field. Similar fluctuations of opposite sign 
are created by electrons, but on much smaller scales that are not resolved in the Figure. 
Accordingly, "striped" zones of positive and negative charge density result, and this elec- 
trostatic analog of a plasma oscillation integrates to produce the E fields in the right panel 
that can accelerate or decelerate electrons and protons. The outcome depends on the shock 
obliquity ©Bni when Obhi ^ 60° , whereas quasi-perpendicular shocks with ©Bni ^ 60° pos- 
sess profiles fairly close to the Gbui = 60° case depicted in Figure 1, since they all have field 
obliquities 9Bn2 ~ 80° - 90° downstream. Note also, that while the gyrational contributions 
are prominent, there is also a diffusive contribution, manifested as an upstream precursor 
modification to p{x) and E. This is particularly marked in the parallel shock ( 0Bni = 0° ) 
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case, where the diffusive scale along the field achieves a maximal component orthogonal 
to the shock plane. This diffusive influence originates in accelerated particles returning to 
the upstream side of the shock {x < 0), enhancing the density there before convecting 
downstream again: protons effect this on larger scales, and so control the precursors seen 
in the Figure (i.e. p{x) > for x < ). Since the fields are established on the scale of a 
proton gyroradius, their magnitude scales as ~ inprg^p = 4:Tienp{mpUixc/ eBi) , so that 
\Ex/Bi\ ~ ATrnpnipUixc/ Bj = M'j^{c/uix) 3> 1 for Alfvenic Mach numbers Ma > 1 • 

The competition between gyrational and diffusive inffuences on electrostatics is a prin- 
cipal conclusion of this paper, defining a dichotomy delineating quasi-parallel and quasi- 
perpendicular shocks. The Monte Carlo technique can accurately trace both influences, 
while comfortably resolving the disparate scales for the e-p shock problem. Since the "lin- 
ear" results illustrated need to be upgraded to account for the E-field's influence on e" 
and p motions, it is presently unclear whether ions can energize electrons overall (the right 
panel of the Figure suggests they may even decelerate them), and how the net work done 
depends on field obliquity. A noticeable feature of the electric field profiles in Fig. 1 is that 
for 0B ^ 60° , these linear field calculations do not establish |E| — > asymptotically as 
1^1 ^ oo 5 as required by net charge neutrality. The next step of this program will be to solve 
the Newton-Lorentz equation of motion dp/dt = g(E-|- v x B/c) to determine both drift and 
accelerative contributions to the charges' motions. These will necessitate a recomputation 
of the E field profiles, and a feedback loop will result, with shock layer currents generating 
magnetic field excursions via Ampere's law, V x B = AnJ/c. This iterative process will 
continue to convergence (establishing |E| — > as \x\ — > oo), with relaxation to equilib- 
rium occuring on the spatial response scale uix/ujp , where ujp — y^Ane^np/nip is the proton 
plasma frequency. Since Uix/{ujprg^p) ~ Uix/{cAiA) ^ 1 ), this response scale is far less than 
a proton gyroradius for typical SNR environmental parameters, and indeed for any strong, 
non-relativistic astrophysical shock. 

The degree of electron energization in the cross shock potential may offer significant 
insights into the well-known electron injection problem at non-relativistic shocks. Electrons 
do not resonantly interact with Alfven waves until they become relativistic. Levinson (1992) 
suggested that e~ interaction with a presumably abundant supply of whistler waves could 
effect pre-injection into diffusive acceleration processes, if electrons could achieve energies in 
excess of around 10 keV to access the whistler resonance branch. The planned self-consistent 
extension of the developments outlined here will help determine whether this channel of 
access to continued acceleration is opened up by shock layer electrostatics. Moreover, crafted 
laboratory plasma experiments might cast light on this aspect of shock layer physics. 
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3. Conclusion 

In this paper, charge density and associated cross-shock electric field spatial profiles are 
presented for different magnetic field obliquities. It was found that in highly oblique and 
perpendicular shocks diffusion plays little role in modulating the field structure, which is 
controlled by the magnetic kinking and compression on the downstream side of the shock. 
In contrast, in quasi-parallel shocks, where the gyrational charge separation is small, diffu- 
sion scales upstream and downstream of the shock dominate the generation of shock-layer 
electric fields. This is an interesting twist, suggesting that observationally, thermal X-ray 
emission could be distinctly different in portions of an SNR rim that establish quasi-parallel 
and quasi-perpendicular shocks. The work discussed here paves the way for self-consistent 
determination of the acceleration/deceleration of electrons and protons, their spatial dis- 
tributions, and the electric fields normal to non-relativistic shocks. This development will 
impact the understanding of electron injection and acceleration in shocks of all obliquities. 
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